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Abstract 

kS \ The correlation functions of the multi-arc complex matrix model are shown 

C^ ' to be universal for any finite number of arcs. The universality classes are 

characterized by the support of the eigenvalue density and are conjectured to 
fall into the same classes as the ones recently found for the hermitian model. 
This is explicitly shown to be true for the case of two arcs, apart from the 
known result for one arc. The basic tool is the iterative solution of the loop 
equation for the complex matrix model with multiple arcs, which provides all 
multi-loop correlators up to an arbitrary genus. Explicit results for genus one 
are given for any number of arcs. The two-arc solution is investigated in detail, 
including the double-scaling limit. In addition universal expressions for the 
string susceptibility are given for both the complex and hermitian model. 
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1 Introduction 

The notion of universality is one of the crucial properties of random ma- 
trix theory. It ensures, that physical quantities of interest such as correlation 
functions or the free energy and its derivatives - the string susceptibility - do 
not depend explicitly on the matrix potential and in this sense are universal. 
The full dependence on the coupling constants can be encoded in few universal 
parameters as the endpoints of the support of the spectral density. 

Universality can be found in different ways when performing the large-A^ 
limit. First, in the so-called microscopic limit the oscillatory behavior of the 
eigenvalue density is studied and found to be universal [ffl] (for recent progress 
see Q and references therein). In this respect particularly the study of the 
unitary ensembles, the complex and hermitian one-matrix model, has drawn 
considerable attention for its application to 3- and 4-dimensional QCD Q. 
There, the merging of the eigenvalue density from two into a single arc provides 
an effective model for the formation of a chiral condensate in terms of random 
matrix theory. 

In the second way of performing the large- A^ limit, the macroscopic limit, 
where the oscillations of the eigenvalue density are smoothed, the splitting of 
the support into several pieces has been much further understood. After the 
seminal works for the one-arc case by Ambj0rn and his collaborators for the 
hermitian and complex matrix model [0, Isl 0], it has been recently shown for 
the hermitian model, that for any number of arcs all correlation functions are 
universal and can be classified completely by the support of the spectral density 

Si- 

The aim of the present paper is to extend the previous results to the com- 
plex matrix model with multiple arcs, generalizing the work of |5|. The same 
achievements will be made here, obtaining again a whole set of universal correla- 
tors characterized by the support of the density. Moreover the new universality 
classes found here are conjectured to coincide precisely with the ones previously 
found for the hermitian model P, 0]. This result was already known for one 
arc ||6[ and is now proven explicitly for two arcs as well. The basic tool therefor 
will be again the iterative solution of the loop equation. Furthermore the string 
susceptibility of both, the complex and the hermitian model with two arcs will 
be calculated explicitly and shown to be universal. 

This opens the possibility to investigate the relationship between multi-arc 
matrix models and integrable hierarchies using loop equation techniques ||9|, |lO| . 
The outcome will have to be compared to earlier results for the two-arc solution 



from orthogonal polynomials |11], which depend heavily on the reliability of the 



ansatz for the recursion coefficients. For the case of two non-symmetric arcs 



this is already very doubtful |12| 



Another earlier attempt towards the universality of the two-arc phase of 
the hermitian matrix model has been made in ||l^. However, the completely 
independent treatment of the two arcs leads to an expression for the correlator of 
densities which fails to fulfill the correct analyticity properties when compared 
to the two-arc eigenvalue density. 



It should be mentioned that macroscopic universahty has been proven for 
orthogonal and symplectic ensembles as well, using variational methods in the 
saddle-point approximation |14]. This result for the correlator of densities has 



been extended to planar multi-loop correlators in the framework of loop equa- 
tions |l^. However, the appearance of odd powers in 1/A^ in the expansion 
make it very difficult to achieve explicit results for higher genera. 

Furthermore renormalization group techniques for matrix models have been 
recently extended to the multi-arc case as well |16|. For the hermitian matrix 



model with a symmetric double- well potential the authors of |16] find a second 
attractive fixed point of the renormalization group transformation, apart from 
the Gaussian one. These results confirm the universality of the planar two-loop 
correlator for two symmetric arcs by completely different means. 

The present paper is organized as follows. After briefly introducing the 
basic definitions in section 2, section 3 deals with the planar solution of the 
loop equations for the multi-arc complex matrix model. Section 4 contains the 
iterative solution of the loop equation, where explicit results are given for genus 
one for any number of arcs. The question of the equivalence with the univer- 
sality classes of the hermitian model is also addressed here. In section 5 the 
two-arc solution is presented in more detail including a proof for the match- 
ing with the corresponding two-arc universality class of the hermitian model. 
Section 6 contains all the possibilities for performing the double-scaling limit 
for two arcs, where the generic case maps to the double-scaled one-arc solution 
ip. Finally, before concluding in the last section the string susceptibility for 
both the complex and the hermitian two-arc model is explicitly shown to be 
universal as well. Throughout the paper the notation for the solution of the 
one-arc complex model M will be followed closely. 

2 Basic Definitions 

The complex one-matrix model is defined by the partition function 

Z[N,{gi}] = e^'^ = y'#td,/.exp(-iVTrF(0V)) , 

oo 

y(<AV) = 5:^(0W, (2.1) 

where the integration is over complex N x N matrices. The generating func- 
tional or one-loop correlator is given by 

W( ^ - 1 V^ (Tr(./'V)") 1 /^, P \ .„ .. 

W(P) = J^l. p2k+i = ]^ \^ p2 _ ^t^ ) • (2.2) 

Introducing the loop insertion operator 



i^*"' --ii^i ''■'> 



the one-loop correlator can be obtained from the free energy, 

W{p) = -^{p)F + i . (2.4) 

dV p 

More generally one gets the n-loop correlator by iteratively applying the loop 
insertion operator to W{p) (or F), 

W( ^ - Arn-2 ^ (Tr(0t0)fei...Tr(0t0)fen)^^^^ 
W{pi,...,pn) = N l^ jj^^ -^j-^^ (2.5) 

= :^iPn)^iPn-i)---^iP2)W{pi), n>2, (2.6) 

where conn refers to the connected part. All the multi-loop correlators and the 
free energy have the same 1/A^ expansion 

oo -. 

Wipu...,Pn) = T.j;^^9iPl^---'Pn), (2.7) 

9=0 
oo -. 

^ = Ej^P.^ (2-8) 

9=0 



SO relation (|2.6| ) holds for each genus separately. In particular eq. (p.4[) may be 
written as 

Wg{p) = ^{p) Fg , 5 > 1 . (2.9) 

3 The Loop Equation 



The loop equation for the multi-arc complex matrix model |17| looks the same 
as for one arc [^ apart from the integration contour C , 

^^ '"^'(^^(u;) = {W{p)f + ^,i-Ap)W{p) ,p^a. (3.1) 



Here C encloses all the square root singularities (cuts) of W{p) along the arcs 
of support a of the eigenvalue density p{y), 



■ [Ji=i{[-X2i-i,-X2i]^[x2i,X2i-i]) seven 

a = ^ , s^ s (3.2) 

( UJi {[-X2i-i, -X2i\ U [x2i,X2i-i])) U [-Xs, Xs] s odd . 

where xi > . . . > x^ > . The integration contour C can thus be depicted in 
the following way 
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Figure 1: C = Uf^^Ci 



As it has been explained already in [18| the eigenvalue density can be intro- 
duced either for real values y or for real positive values y^, y € R, where y^ is 
the eigenvalue of 4)^4). In the first convention, which is chosen throughout the 
calculations, p{y) is a smooth, symmetric function on R, whereas in the latter 
it becomes singular when having support at the origin |jl^. In eq. ( |3.2| ) this is 
just the case when s is odd. The corresponding density on R+ is obtained by 
transforming /9(y) -^ p{y)/2y = p{y'^). 

Now there is a simple way to see why the case of an even s is actually 
equivalent to the hermitian model with s/2 cuts on R+Q. When identifying 



-W(p) 
P 



Xi > , i = 1, 



,■5 , 



(3.3) 



one exactly arrives at the loop equation of the hermitian matrix model (H) with 
s/2 cuts Q, plus an additional term ~ 1/p. This term must vanish since W^ {p) 
has to be regular at p = for s even, as the origin is not included in the support 
a thenQ. In |18| the same observation had been made on the level of orthogonal 
polynomials for the case of s = 2, identifying the respective string equations of 
the complex and hermitian model. Consequently in all the following only the 
case of s odd will be further considered. It will be called {s + l)/2-arc solution, 
counting the number of arcs on R_|_. This way of counting is chosen to name 
the same universality class in the hermitian and complex model by the same 
number of arcs. 

Solving the loop equation (|3.l| ) in the planar limit the 1/A^'^-term can be 
omitted. The result is given by [|l^] 
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(3.4) 



^ On the level of the eigenvalue model in |18| this is intuitively clear. 

^ This fact had been missed i n_J17| , giving a simple explanation why the results presented 
there for s = 2 are equivalent to |6|^ 



which can be rewritten as 



W,{p) = \(v'ip)-Mip)Jl[ly-xj)) . (3.5) 



4 
Here M(p) is an analytic function given by 



,,, , f duj 2tjy'(tj) 1 ,„^, 

M{p) = i IT-. 2 i I ■ 3.6 



The signs of the square root in the complex plain in eq. ( |3.5| ) are chosen to be 
{^pY at iboo. From eq. ( ^I^ ) the eigenvalue density p{y) can be immediately 
read off: 

p{y) = lira--—{Wo{y-ie)-Wo{y + ie)) ,y£a, 



= i 1^(^)1 v-ni>^^- (3-7) 

The parameters xi, i = 1, . . . , s, given as functions of the coupling constants 
Qi, i € N+, are determined by the boundary conditions on Wo{p). From the 
asymptotic behavior W{p) -^ ^ eq. (|2.2|) it follows 

1 f duj I , , 1 



^'- = 2 fc 4^" ^ (") M. , , ., , / = 1, 3, . . . , . . (3.8) 

In analogy to the hermitian case Q the remaining ^^ equations for the Xi are 
obtained by a criterion of stability |p^P| 

= / dyM{y)J\l iy^-x^^) , fc = 1, . . . . (3.9) 

For s > 3 these equations lead to the explicit appearance of elliptic integrals 
in the higher genus correlators as well as in the planar n-loop correlators for 
re > 2. 



4 The iterative solution 

After having solved the loop equation ( |3.1| ) in leading order, an iterative 
scheme can be built up for calculating higher order contributions. It is based 
on the planar solution ( |3.4D along the same lines as in |^], which is treating the 
one-cut situation. Inserting the genus expansion of the one-loop correlator eq. 
( p.7| ) into the loop eq. ( |3.1| ) and comparing order by order in genus, the higher 
order contributions are determined by 

{it - 2W^{p))Wg{p) = Y. Wg:{p)Wg-g,{p) + —{p)Wg^^{p) , g>l. (4.l) 

g/=l 



Because of symmetry the stability criterion has to be imposed only between half of all 
arcs. 



Here A^ is a linear integral operator, given by 

JC iTTl p^ — UJ^ 

Once the result for Wq (p) is known one can determine Wg (p) for g > 1 iteratively 
from contributions of lower genus on the r.h.s. of eq. ( |4.1| ), provided that the 
operator {IC — 2Wo{p)) can be uniquely inverted. However, it can be shown 
from the boundary conditions ( |3.8| ), that this operator possesses the following 
zero modes 

KeT{IC-2Wo{p)) = Span{,/.W(p),pV^o)(^)^___^p.-yo)(^)|^ 

^(0)(p) ^ 1 (4.3) 

^ULiip'-^l) 



Because of the asymptotic of Wg{p) ~ ©(-j) for g > 1, which can be seen from 

the definition ( |2.2| ), the last zero mode ~ - cannot contribute to Wg{p). In the 
case of s = 1 this argument uniquely fixes the inversion of eq. ( |4.1| ) ||5|. Here, 
the condition ( |2.9D will be used in addition to achieve a unique solution Wg{p) 
for any s. Therefore a set of basis functions is introduced in the following way: 

{JC-2Wo{p))x^\p) = -^^^-, i = l,...,s , nGN+, 

\P ^i ) 

{JC - 2Wo{p)) V^")(P) = 4^ , n G N+ . (4.4) 

In addition a change of variables from the coupling constants gi to the following 
moments is performed: 



/. 



(fc) 



duj ll>V'(lo) ,rm/ X • -, , TVT 



C 47ri {uj"^ — xJY 



Apart from the endpoints of the cuts x, these moments will encode the univer- 
sality of the higher genus correlators. Namely given that the r.h.s. of eq. ( [4.1| ) 
is a fractional rational function oip having poles at the xi only, Wg{p) will have 
the following structure: 

3g— 1 s g 

Wgip) = EE^t^^^^P) + E<^V'("^(P)- (4-6) 

ra=l i=l 71=1 

The coefficients of the basis functions A^'^J and Dg' are complicated functions 
of the endpoints Xi and of the moments eq. ( |4.5| ). As the order of the highest 
pole in Wg{p) is not changed compared to the one-cut case ||5|, Wg{p) will 
depend on at most s{3g — 1) moments I^^ and g moments M'"^ This follows 
from the same arguments as in 0]. The unique set of basis functions can be 



proven to be (see appendix^): 



xf^(p) = -j^ 



bt\p) 



n-l N 



ip(^\p) 



1 



M(o) \ 



[7(")(p) 



:^yT-part 



1 



(p2 — xfY' 



k=l / 

n-l 

fc=l 



1, . . . ,s , n G N_| 



^^"^(^) = 4r'^^°^(p)' 



n G Nj 



(4.7) 



The restriction to the part which can be written as a derivative of the loop inser- 
tion operator -ijiji) ensures that the whole p-dependence of the basis functions 
and thus of Wg{j)) in eq. ( |4.6| ) can be absorbed into such derivatives in or- 
der to fulfill eq. (|2.9D. At the same time this requirement uniquely fixes the 



ambiguities coming from the zero modes in eq. (4.3). 

The explicit construction of the basis functions can be found in appendix 
A, the first of which reading 

1 dx? 



xf^(p) 



X?^) 



V^(i)(p) 



X? dV 



(P), 



-2 dL 



(1) 



3xf/f ) AV 

2 (iM(0) 
~M(o) dV 



ip) 



= l,...,2s , 
2 dx] 



3xf dV 



W-s^S 



'dx'i 



dx'i 



2 



jC^f^x-i^w 



" 1 r/r^ 



■ 1 xf dV 



(4.^ 



The one-loop correlator of genus one Wi{p) can now be calculated from eq. 
( [4.1| ) for 5 = 1, which reads 



(^-2Wo(p))Wi(p) = -^{p)W^{p) 



(4.9) 



Therefore the planar two-loop correlator Wq{p,p) on the r.h.s. has to be deter- 
mined, which can be achieved from applying -^{p) to eq. (|3.4| ) in the following 
form 
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/dxf 
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dxf 
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OO ^^(fc) 
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dy 
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fc=0 



dV 



(4.10) 



Using the results from the appendix A eq. ( |A.7| ) one obtains 



Wo{p,p) 



1 '^ 

-E 



x: 



1 "^ 

-E 



x: 



1 



16 ~^ {jp^ — xfy 16 f-^^ xf — x^ I p2 _ xj 

s-3 



p2 — X^ 



+ 



1 s A 1 1 A v^ a,^i xf 



16p^ 



16 ^p2 



+ 7EE 



4 ^ ;^ p2 _ ^2 ' 



(4.11) 



where the ai^i are given by the solution of the hnear set of equations ( [A.^ ). 
They depend only on the endpoints of the cuts Xi and on elliptic functions of 
them. The result for Wi{p) now follows from inverting eq. ( |4.9| ) by inserting 
the basis functions eq. ( |4.8D into Wq{p,p): 



Wiip) = ^t^h?\p) - ^tzS:;^{x^\p)-xf\p) 






s-3 
S 2 



+ l^V^^'^(^)-^Exf^(p) + ^EE«M^?xf^(p)- (4.12) 

The calculation of higher genera Wg{p) as well as the explicit integration to 
get the free energy Fg becomes very technical as one has to differentiate or to 



integrate the set of elliptic functions given in eq. ( A. 6 ). Still, following the 
arguments from |^ it becomes clear, that the general structure for Wg{p) will 
be precisely the one given in eq. ( |4.6| ) for an arbitrary number of cuts. 

In the case s = 3 the elliptic integrals occurring are given in terms of the well 
known complete elliptic integrals of first, second and third kind (see eq. ( |B.2| )) 
making more explicit results accessible. In the generic double-scaling limit the 
solution will be mapped to the double-scaled one-cut solution [y]. This relation 
is probably true for all multi-cut solutions. In [^ higher genus results have been 
obtained explicitly up to and including 5 = 4. 

Apart from the general result for genus one for any s also the planar two- 
loop correlator at different arguments can be calculated explicitly, using again 
the results from appendix A. Performing a similar calculation as for Wo{p,p) it 
follows 

^OW,P) 4</>(0)(p)l^(g2-p2)2 + g2_p2l^ 1 + Z.p2_^2 

s-3 

Thus WQ{p,q) and because of eqs. ( ^ ) and ( [4.10| ) also all the higher planar 
multi-loop correlators are universal in the sense that they only depend on the 
endpoints of the cuts Xj. The same statement can be made for the density cor- 
relators, as they can be obtained from the multi-loop correlators and vice versa 
(see e.g. Q). Consequently eq. ( 4.13 ) gives a classification of all universality 



classes of the multi-arc complex matrix model in terms of the support of the 
eigenvalue density. 

A similar classification has been found in previous works for the hermitian 
model 0, P] , so a natural question to ask is whether one can identify the classes 
here and there, which have the same number of intervals as support. In order to 
do so one has to map again the correlation functions from the hermitian model 
Q to the corresponding quantities in the complex model as in eq. (pl^), 

\wo{p,p) ^ Wi'{p,p) _ , . (4.14) 



p2 



p=p^ 



In addition the smallest endpoint X2s of the cuts of the herniitian model {H) has 
to be set to zero, as now the complex model with odd s is treated^. While for 
the one-cut case it is known already, that the universality classes for the planar 



two- loop correlators match y, in the general case Wo{p,p) from eq. ( [4.1l| ) has 
to be compared with Q 

16 ^ [p- XiY 16 f-^^ Xi - Xj \p-Xi p- XjJ 

2s-4 
1 o- 1 2s— 1 -I -I 2s— 1 2 - -J 

+ —- ^^"^ V + i V T "^-'^^ + ^^(4 15) 

I6p 16 ^p-x, 4^;e^p-xi 4p ^ ■ ^ 

for p = p^ , Xi = X?, i = 1, . . . , 2s — 1 and X2s —^ 0. When identifying s = 2s — l, 
which is the corresponding number of endpoints, the two correlators eq. ( [4.11| ) 
and eq. ( [4.15| ) are identical apart from the last extra term in eq. ( [4.15| ) which 
is non vanishing in the limit X2s — > 0. However, one has to bare in mind that 
the Ui^i are solutions to a set of 2s {s — 1) linear equations (see eq. (5.16) in |^]), 
which are different from the s{s — l)/2 eqs. ( [A.8| ) determining the a-i^i. A proof 
for the equivalence of the two sets of universality classes would thus need some 
nontrivial identities for the elliptic integrals occurring and is therefore left as 
a conjecture. Still, for the case of s = 3 this proof can performed explicitly as 
it will be shown in the next section. Summing up the results obtained so far 
the complex model with both s even and odd contains all universality classes 
found for the unitary ensembles, due to the fact that the hermitian model is 
contained in the complex one for the case of s even. 

5 The two-cut solution s = 3 

As it has been mentioned already for s = 3 further results become accessible 
as the elliptic integrals originating from the boundary conditions can be related 
to the complete elliptic integrals of the first, second and third kind. Explicit 

expressions for the zero mode coefficients Oj of the quantities -^ (p) , which are 
responsible for the appearance of the elliptic integrals, are given in appendix B 



eq. 

In the previous section the set of universality classes for all multi-arc so- 
lutions for the complex model has been given and an equivalence to the ones 
recently found in the hermitian model 0, ^ was conjectured. In the special case 
considered this equivalence will be shown to hold. For s = 3 the general results 



obtained for Wo{p,p) eq. ( 4.11 ) and Wo{p, q) eq. ( 4.13 ) take the following form 



using the results for the ai from appendix B, 



*In matching the correlators of the two models the density for nonnegative eigenvalues of 
the complex model is considered. 



1 A xf 1 1 A xf + x'^i 



Wo{p,p) ifiZ^ ^^2 _^2^2 + lfi^2 irZ^ (^2 _^2^|-^2 



16 ^ (p^ — x?)2 ' 16p2 16 ^ (p2 — x?)(p2 — x?"! 



1 1 3(xi + xi)-2p2 1 xi(x2-xi) i^(A:) 



Wo{p,q) 



I6p^-xl 16(p2-xi)(p2-a;2) 4nf=i(p2-x2)K(A:) 

(5.1) 



xl){q'^ - xl){q^ - xl) 




X|)(p2-x2)(g2_a;2) 



2^/'^2 



x|)(p2 -xl){q 



xDiq"^ - xDiq"^ 

E{k) _ 1 

K{k) ~ 2(p2-g2 



+ 7 ^^===^^=^i^ - ^TZ^"^^ ■ (5-2) 



The consistency of the two results can be checked exphcitly. Namely Wo{p,q) 
is symmetric, it has the same analyticity properties as Wo{p), it asym,ptotically 
reaches ~ 0{^), it is regular at equal arguments and coincides with Wo{p,p) 
and it finally has the correct discontinuity — , i^i-^i along the cuts. The con- 
jectured equivalence of the universality class of the hermitian model with two 
cuts at [x4,X3] U [x2,2;i] is now evident when comparing with the results given 
in Q (eqs. (6.6) and (6.5)). Setting X4 = and identifying Xj = xf it holds 

—Woip,p) = W(f(p,p) p^p2 

-Wo{p,q) ^ T^o''fe9)|p=p2.,=,2 . (5.3) 

■^^ XA^X., X/1 — 

So despite the fact that in the complex model the eigenvalue density p{y'^) 
becomes singular at the origin whereas the density of the hermitian model 
stays smooth, both models belong to the same universality class for one and 
two cuts. Nevertheless, for higher genera the two models will give in general 
different results, as it can be seen for example from comparing the free energy 
obtained below with the one in g. It is only in the double-scaling limit, that 
all the models mentioned will become completely equivalent (see sect. |6|). 



The correlator of densities given in |13|, which can be related to the planar 



two- loop correlator WQ{p,q) (see e.g. Q), does not match with the above two- 



arc universality class for the following reasons. First of all the authors of |13] do 
not impose any boundary condition as in eq. (|3.9| ) and are therefore left with a 
parameter dependent expression. Furthermore their independent treatment of 
the two arcs, adding up two single arc correlators with different endpoints, does 
not fulfill the correct analyticity properties compared to the two-arc eigenvalue 
density in p. 

Turning back to the results for genus g = 1 from the previous section also 
the free energy Fi can now be determined. This is due to the fact that the 



10 



derivatives of the complete elliptic integrals with respect to their moduli are 
known and close among themselves (see e.g. [PQ]). From eq. ( 4.12| ) the one- 
loop correlator of genus one for s = 3 takes the following form, 

1 ^ 1 fdxf .xldxl . 1 1 rfMW 



16 ti" ^2 - xH dy ^'^ x^: dV ^^'J 8 M(o) dy 

1 .^-^ 1 dxf . . I y-^ tti dxf . . . ^. 

2 = 1 t ?.= 1 *■ 

Using the results for the a^ eq. (|B.5|) and standard formula for the complete 



elliptic integrals |2C] Fi can be integrated according to eq. (2.9) 



3 3 

i=l i=l 

1 1 

-Qj2^'^\xl-x]\ + -{ln\xl-xl\+ln\xl\) . (5.5) 

i<j 

Apart from the factor of the M'-'^-'-term all the other terms coincide with the free 
energy of genus one of the hermitian model |^] at X4 = 0. The same observation 
can be made when comparing the free energies of the one-cut solutions. There, 
at higher genera the difference between the hermitian and complex model be- 
comes more evident. The iterative procedure of solving the loop equation is 
purely algebraic and in the case of s = 3 even accessible to computer algebra, 
grace to the knowledge of the elliptic integrals. In the next section higher gen- 
era will be obtained by other means, namely by comparing to the results for 
the one-cut solution in the double-scaling limit j^. 

6 The double-scaling limit for s = 3 

In general two different situations may occur in the double-scaling limit 
(d.s.l.). In the generic case, when the limit is performed at any of the endpoints 
ibxj, the contributions coming from the multi-cut structure will drop out and 
the correlators will precisely match with the continuum version of the one- 
cut case [^]. A different situation is given, when the d.s.l. is performed at 
simultaneously shrinking or/and merging arcs. The correlators will differ from 
the one-arc result, as it has been found similarly in the d.s.l. of the hermitian 
two-cut model ^. 

6.1 The scaling limit at Xj 

Choosing any endpoint, say Xj, the m-th multi-critical point is defined by ?7i — 1 



extra zeros of the eigenvalue density eq. (3^) accumulating at +Xj and 



11 



By tuning the coupling constants in order to reach this point the argument of 
the correlators p and the ±Xj will scale like 

p^ = {x'-jf + avr , 



X 



„2 (^c\2 



3 



(x5)2 - ak^ , (6.1) 



p=Xi 



with a ^^ 0, whereas the xf,- do not scale. Since l\ = ^ o^^iy ^ 2{k-i) 
the j-th moments will scale like 

if^ ~ a™'^' , (6.2) 

whereas the other moments l\J- and M^") do not scale either. In the limit 



eq. (|6.1| ) the elliptic integrals K{k) and E{k) and consequently the zero mode 



contributions aj will stay subdominant in eqs. (A_/?), ( [A.3| ) and ( [A.4| ), 
^{P) = -^^f\p) (d.s.l.), 

3 



(6.3) 



^y{p) = {k+-)\lY^^>^{p)-x'^<t^';'^^>{p)\ (d.s.l.). (6.4) 

Only the basis functions Xy {p) ~ a"'"""^^ will give dominant contributions 
when inverting the loop equation. Moreover the loop insertion operator itself 
will simplify due to eq. (6^) as well as Wq{p,p) in eq. (|4.11| ) , the starting 
point of the iteration, 



^3> 
^2 a oo ^r(*=) 



16 {p^ — X 



d dxi d ^ dif' d 

3 fe=l '-'^j 

The eqs. ( 16.3|) - (|6.5D precisely map to the corresponding quantities of the 
double-scaled one-cut solution, which are given in the appendix A of Q. Con- 
sequently both solutions are equivalent in any order of genus, as they start the 
iteration from the same quantities. In fact more was proven in Q, namely that 
in the d.s.l. the complex and the hermitian matrix model become completely 
equivalent. From the above results this equivalence now extends to the complex 
model for s = 3. The same matching has been shown for the two-cut hermi- 
tian model [p. This equivalence has to be understood in the symbolical sense 
that the basic objects as moments and endpoints of support appear exactly in 
the same way. However, when expanding these objects in terms of coupling 
constants gi, the expansion explicitly depends on the number of cuts. 

6.2 The scaling limit at merging and shrinking arcs 

In the case where the d.s.l. is considered at an endpoint of an arc, which is 
simultaneously shrinking to zero, merging with another arc, or both together, 
the following situations are possible: 
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ij X2 -^ Xs merging 

ii) X3 — > shrinking 
iii) X2 — > xi shrinking 
iv) X2,X3 ^ merging + shrinking 

v) xi,X2 -^ X3 merging + shrinking 



-Xi -X2 -Xs Xs X2 Xi 



Figure 2: The support a of the eigenvalue density p{y) for s = 3 

In ah five cases the resulting continuum limit will differ from the previous 
one called generic, and thus from the d.s.l. of the one-cut solution. Only in the 



last case v) the contributions containing elliptic integrals in eqs. ( A. 7 ), ( [A.3| ) 
and ( |A.4| ) will stay dominant. In all five cases the iterative solution of the loop 
equation with not simplify much any longer. The case i) will be treated in more 
detail as an example, whereas in the remaining cases only the scaling behavior 
of the zero mode coefficients Oj of the quantities -^^(p) will be given. 

i) The transition of merging cuts can be parameterized by 

2 2 I 

X2 = x^ + av , 

a;| = x^ — a/i , z^, /i > , 

p = x^ + aiT . (6-6) 

In this limit the modulus k of the elliptic integrals vanishes, 

k' = a{u + ^,) f + 0{a') , (6.7) 

and consequently K{k) and E{k) can be expanded in terms of /c^, 

^(^) - 1 _ 1^2 + Oia^) . (6.8) 



K{k) 2 

The coefficients Oj will stay finite in the limit 



"1 = ^2_^2 + ^(«) ' 

"^ = I - 2{xl- xl) +^^")' 

03 = aa + 0{a) , (6.9) 
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leaving the zero mode contributions aicl)^^'{p) ~ a ^ subdominant com- 
pared to (t>j=2 3 ~ '^"^ ill 6q. (|B.3|) . Due to the scaling of both sets of 
moments /• Jg 3 their derivative now reads 

(6.10) 

and the same for the indices 2 and 3 interchanged. The starting point for 
the iteration Wo{p,p) becomes 

W ( \ - L\- xf 1 (xj + xj) / 1 1__\ 

Wo[P,P) - ^g ^2_.^ ^^2 _ ^2^ ^g ^^2 _ ^2) [p2 _^2 p2_^2) id.S.i.j 

(6.11) 
from which Wi{p) can be easily obtained. Clearly the result differs from 
the generic d.s.l. in the previous subsection. 

ii) The d.s.l. at the inner arc shrinking to zero can be parameterized as 



a 



p^ = air . (6.12) 

As the modulus k'^ reaches unity in this limit the elliptic integrals can be 
expanded in terms of the complementary modulus k''^ 



r2 ^2 



/2 _ 1 7.2 _ „ ^1 ^2 



1-^2 



The ai will stay finite or vanish, 

ai = as = 1 + 0{- — ) , 

ma 

as ~ 0{-^) . (6.14) 

ma 

Still, the continuum limit will be changed compared to the generic d.s.l. 
as now the moments M' •* will also scale. 



iii) The shrinking of the outer arcs is parameterized by 

2 2 , 

Xi = x^ + ai' , 



^2 — 3;(, — a/U , fi,iy > , 

P = X^ + UTT . (6.15) 
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Expanding again in the complementary modulus k'"^ leads to 

for the elliptic integrals and hence to the following behavior for the zero 
mode coefficients 

ai,a2 ~ C 



a In a' 



Although they become singular in the limit ( |6.15| ) they still remain sub- 
dominant, leading again to eq. ( |6^ ) for j = 1,2. The starting point 
for the iteration becomes like eqs. ( |6.10 ) and (6.11), with the index 3 
replaced by 1. 

iv) The merging of the arcs while the inner arc is shrinking to zero is para- 
meterized by 

X2 = av , 

X3 = a/i , i^ > // > , 

p^ = avr . (6.18) 

In this limit the modulus reaches 

fe^ = 1 - ^ + 0{a) , (6.19) 

V 

which is not necessarily small and may take any value in (0,1). The a, 
become 

ai = 1 + 0{a) , 

"3 = ^— f-/" + ^^^ttt) + C'(a) . (6.20) 



Consequently they do not contribute in the d.s.l.. 

v) In the case of merging arcs while the outer arcs are shrinking to zero the 
ai will be sufficiently enhanced. The d.s.l. reads 



2 


= x^ + au , 




xl 


= xl + afj, , 1/ > fi > , 




4 


- x^ 




dull 


= x1 + a-K . 
IS becomes 


(6.21) 




k" = ^ + 0{a) , 


(6.22) 
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which can reach any value in (0,1). The zero mode coefficients now read 



a2 = a 



as = a 



In this case all Qj become sufficiently singular to contribute in the -^{p) 
in the d.s.l.. The iterative procedure will therefore not simplify and in 
addition the elliptic integrals remain explicitly present in the continuum 
limit. This interesting new continuum behavior has also been found in 
the hermitian model [01. 



7 The string susceptibility 

In this section a closed and universal expression will be given for the string 
susceptibility of both the complex and hermitian matrix model with two cuts. 
It reflects once more the fact that the two-point correlators of both models are 
universal as well. As the calculation for the complex model is much simpler 
it is presented in some detail. For the hermitian model only the final result 
will be given, which can be obtained along the same lines when starting from 
H. Both susceptibilities are shown to belong to the same class of universality. 
The knowledge of the string susceptibility may either serve to study the critical 
behavior of the model, as for example in |2l| for the 0(n)-model or in [^2[ for 
the Penner model. On the other hand when one tries to relate matrix models to 
integrable hierarchies using loop equation techniques |^ 10 1 the susceptibility 
enters directly the corresponding Gelfand-Dikii ansatz. 

In order to introduce an overall coupling constant in front of the potential 
the following replacement is performed, 

V{p) - ^Vip) . (7.1) 

The parameter T can be thought of as a temperature or the cosmological con- 
stant. The string susceptibility U{T) is then obtained from the planar free 
energy Fq, 

U(T) = -|j(r=F„) 

~ (Tc - T)-^''"- . (7.2) 

Here it is supposed, that U{T) shows some critical behavior at T = Tg. In order 
to exploit the knowledge of the quantities calculated so far as Wo{p), the two 
following identities will be used, which can be shown easily from the definitions, 

rr2dFo _ ^ I duj 



4T -/c4„:''"°(-'^(-'- 

r^ = -4 (t^'-^) =^A2i ^TW,i.)n.) ] . (7.3) 

dT dT^ \ dT J dT^ \ Jc^TTi "^ ^ ^ ; ; v ; 
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The final result obtained will thus be ^. Taking Wo{p) from eq. (3.4) together 



with the replacement (7J.) the first differentiation of the integrand yields 



dx'^ 

In order to proceed the quantities -^ have to be determined, which can be done 

exactly along the same lines as the calculation of the -^ (p) . Here s = 3 has to 
be specified to obtain explicit results. Differentiating the boundary equations 



( p^ ) and (3^) with respect to T leads to the following set of equations 



^ ^ dT 
1=1 



1=1 



T^Tll'^^K,- (7-5) 

4 = 1 



The solution is obtained with the help of appendix B eq. ( [B.2| ), where the K, 
are given, 



dT U,M^l-x])V'' ""' K{k) 



Tir-^ = rr f^2 .A ^'-4^7^) , i = 1,2,3. (7.6) 



Inserting this into eq. ( |7.4| ) for s = 3 one arrives after some calculation at 

A(™„M)^*.».M(.^-.?ffi||)), ,„) 

and consequently 



As the differentiation of the complete elliptic integrals with respect to their 
moduli is known (e.g. in [^), it is useful to rewrite ^ in the following way 

d d -A^ dxj d dk"^ d do? d 

1=1 * 

Using this form on the r.h.s. of eq. ( f7.^ ) one arrives again after some tedious 
calculation at the following result for the integrand 



A(^<^(0)(^)(^2 .2n(a,fc) 

3 j^2 / ^2 



"^'^W 



igt(-|^)(^<^i"H + -*'"'H). (7.10) 



In order to get a universal expression for -^ in eq. (^j^) one would like to 
perform a partial integration to arrive at an integral containing V'{uj). Next the 
boundary conditions eqs. ( |3.8[ ) and ( |3.9| ) would allow to eliminate the explicit 
dependence on the potential. Indeed the wanted transformation can be done 
due to the following observation. The last parenthesis in eq. ( |7.10| ) containing 

the w-dependence is precisely given by II -^{oj), which can be written as a 
total derivative with respect to u). This becomes evident when rewriting the 
linear set of equations ( |B.1| ) , which determine these quantitiesF], 

1=1 " 

3 ,_,^^2 



E-?/f^^(p) + |(p¥°)(p)), 



1=1 

1=1 " 

A{P) = / dy-^ 21M7T ('^•1^) 

The last equation can be shown to hold using eq. ( |A.6| ). In the same way the 
boundary condition (^^ is shown to be equivalent to 



= (( ^V'{uj)2uj(t)^^\uj)A{Lo) . (7.12) 

Jc 47^^ 

From eq. ( [7.11| ) one obtains 

(7.13) 
where 

E = {xl- xl)Ki + {xl - xl)K2 + {xl - xl)K3 . (7.14) 

dx'^ 

The -^(p) for i = 2, 3 are obtained by a cyclic permutation of the indices. 
When inserting eq. ( [7.13 ) into eq. ( 7.10 ), integrating by parts in eq. ( |7.S| ) and 

iisincr tlip Kminrlarv rnnrlitirms (li fi\ anri (7 ^'?.\ nnlv tVip tprms r~~^ iT'di^^) (m\ 



using the boundary conditions (|3.8| ) and ( 7.12 ) only the terms ^ p (p^ '{p) will 



survive in eq. ( |7.13| ). The result for the string susceptibility finally reads 
dU _ ^ dxj 2 / , _ 2 n(a,fe) y 

It is a remarkable universal quantity, as it shows no explicit potential depen- 
dence. Compared to the one-arc susceptibility 



d\Js=i 2 dx' 



dT j;2 ^rj, 

dx^ 4 



dT TI(^) ' 



(7.16) 



The same can be achieved for the general set of equations ( |A.5| ) and (|A 
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which can be obtained as a simple exercise, eq. (|7.15| ) contains much more 
structure and thus offers more possibihties for showing critical behavior. 

The string susceptibility for the hermitian matrix model with two cuts can 
be calculated exactly along the same lines presented in this section, taking ||8| 
as a starting point. It is given by 

where 

(7.18) 
To compare with the result for the complex model eq. ( 7.15| ) X4 is set to 
zero, as being done already in section Q. It can be seen now explicitly, that 
both string susceptibilities belong to the same universality class. This is not so 
astonishing, since the planar two-loop correlators of the complex and hermitian 
two-cut model have been identified already in section ^. For completeness the 
one-arc susceptibility is also given for the hermitian model, 

dU^^^ _ ^ dxi 2 



dT .f^2 ^^ Ylj^ii^i - X. 



]j 



-(1) dxi 4 



TD''—^ = -, i = l,2. (7.19) 

8 Conclusions 

In this paper the program of classifying and calculating correlation func- 
tions of the unitary ensembles, the hermitian and complex one-matrix model, 
is completed. For any number of arcs of the support of the eigenvalue density 
the correlation functions are found to be universal for both models, thus pro- 
viding a complete classification of all possible correlators to occur. The higher 
order corrections in l/N"^ to the correlators are determined by an iterative so- 
lution of the loop equation, where for genus one explicit results are obtained for 
any number of arcs. In general, in the double-scaling limit the one- and two-arc 
solution of the complex and hermitian model become equivalent all together. 
However, when the scaling limit is performed at merging or /and shrinking arcs, 
a new continuum behavior has been found here in the complex model as well 
as in the hermitian model Q- 

An interesting question is to find applications for the multi-arc solutions, 
which were calculated here in the macroscopic large- A^ limit. For example in 
solid states physics such situations with a Hamiltonian having a band structure 
for its eigenvalues typically occur. 

It seems to be straightforward to analyze the critical behavior of the two- 
arc solutions now directly from the given string susceptibilities. This could 
be done in a similar way as the investigation of all possible transitions in the 
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double-scaling limit in section ^. Another point left for future work is the 
relation between the multi-arc matrix models and integrable hierarchies. The 
use of loop equation techniques [^, IC] requires the knowledge of the string 



susceptibility as it enters directly the corresponding Gelfand-Dikii ansatz. In 
pA[| the mKdV and a subset of the nonlinear Schrodinger (NLS) hierarchy has 
been proposed for the symmetric and non-symmetric two-arc solution of the 
hermitian matrix model. However, in the latter case the ansatz made there for 
the orthogonal polynomials is very doubtful |12], so it remains to be checked 



whether the same results can be obtained from the loop equations. 
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A Derivation of the basis 

It can be easily shown by induction, that the functions given in eq. ( [1.7] ) 
without the restriction to the ^-part form a possible set of basis functions, 
satisfying eq. ( ^^ ) up to the addition of allowed zero modes of eq. ( ^ ) : 

, i = l,...,2s , n e N+ , 




M(0) 



[7(«)(p)_^^(fc)(p)M("-^) , nGN+. (A.l) 

fc=i / 



In order to absorb the p-dependence of these functions into derivatives of the 
loop insertion operator, the ^(p)-part of (pf" {p) and i7(")(p) has to be calcu- 
lated. The basis functions are then fixed uniquely by subtracting the remaining 
part, which will turn out to consist only of zero modes. Therefore the -^- 

derivative of the basic variables l) , M^^> and xi has to be determined. This 



can be achieved by rewriting the loop insertion operator as in eq. ( 4.10 ) and 



applying it to the definitions (4.5). Using the identity 
the results read 



^AV)V\-) = IZr^Tvl ' (A-2) 
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dV 



(P) 



2k + I 



^^'^'\P)- '^t4T2{-M'Hp)-t'§ip) 



'^t^i^i 



■l±i§ip)t::^ - t:j(^^'Hp)]-iAA) 



i=l 



1=0 Xi 



1=1 Xi 



dxj 



The -^ (p) can be obtained from solving a linear set of equations derived from 
the boundary conditions. Applying ^(p) in the form of eq. (4.10) to the first 



set of boundary equations (3^) leads to 

= J2 (xflP^ip) - ^.V-Vf ^(P)) +{k-s) p'-'<P^'Hp), k = l,3,...,s 

(A.5) 



i=l 



For the second set of boundary equations (^) one get^ 

X2j 







2dp 



dy 



X2j + l 



2p (j)(^\p) 

p2_y2^(0)(y) 



E4^^#(^)^m 



i=l 



dV 



T.[-l^?'ip) -4^^^(p))i^M, 



J = 1, . . . 



s-1 



K 



X2j 

i,j = I dy 

X2J+1 



Making the ansatz 

„2 



dV 



m=i{y'-4) 

(y2 _ xj) 



s-3 
2 



(A.6) 



^i^(p^ = ^"^i^p) + E«.^^'V(°)(P) 



i = l,...,2s 



(A.7) 



1=0 



where the sum is over the allowed zero modes in eq. (4.3), one arrives at the 
following set of equations, 



3-3 

2 



,1 xf /' - E ^T'^ ^0 + (2A; + 1 - .) p2^ fc = 0, . . . , ^ 

i=\ 1=0 1=0 



E ( E «i>' 



s-3 
S 2 



= E E «».' ^ij p^^ ' ■?■ = !' 

i=l 1=0 



(A. 



Comparing powers of p leads to precisely s(s — 1)/2 linear equations determining 
the ai^i, i = 1, . . . , s, / = 0, . . . , ^^, completely as functions of Xi and Kij, 
which are elliptic functions of the Xj. The Oi^i are thus uniquely determined 
by the support of the eigenvalue density a. In appendix B they are calculated 
explicitly for the example s = 3. 

After having determined the -^{p) the -jy-part of 4>\ (p) and Q^'''{p) for 



the basis can now be obtained as follows. Resolving eqs. ( [A.3| ) and ( |A.4| ) for 

"The argument p has to be excluded to be between the cuts during the calculation and 
analytically continued in the end (compare to |8|, appendix A). 
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dk+l) 



(p) and Q^^^'{p), it becomes evident together with eq. (|A.7|), that the 



dv^iP) oi^ -dv^iP) IS just a 



part which is not expressed in terms of -j^ip) 
Unear combination of allowed zero modes. These terms can now be subtracted 
in eq. ( [A.l ) in order to achieve the unique basis eq. ( [4.7D . The firsts terms 
needed for the results in eq. (|4.8| ) read 



/,(!) 



/.(2) 



ip) 



ip) 



-rrr-part 



-part 



Q^^\p) 



-jTT-part 



r(i) 



dxi 



xf dV 



(P), 



l,...,s 



iP dx^ 



(1) 



xf dV 



Hp) 



J^dl\_ 



■ip) 



2I^dxl 
3xf dV 



(P) 



r(i) 



4 = 1 * 



dV 



(A.9) 



da;. 



B The ^{p) for s = 3 



dV 
In this appendix the equations determining the zero mode coefficients Oj^/ 

of the quantities -^{p) are explicitly solved for s = 3. To start with the set of 
equations (|A.5|) and (|A.6|) for s = 3 read 







j:[ii'^^iP)-xu^\p)^-2cb('Hp), 






/,(0) 



(p) E ^? ' 



i=l 



a;2 



iCi 



dy 



X3 



y2 - xf 



(B.l) 



After substituting y^ —f y in the elliptic integrals Ki they can be found in [ |20| , 
where they are expressed by the complete elliptic integrals of first, second and 
third kind: 

Ki = x[xlxlK{k) + xlixl-xl)E{k)+xl{xl-xl-x^s)n{a,k)) 

K2 = x({2xl-xl)xlK{k)+xl{xl-xl)E{k)+xl{-xl + xl-xl)Ii{a,k)) 

K3 = x[{2xi-xl)xlK{k)+xl{xl-xl)E{k)+xl{-xl-xl + xl)n{a,k)) 



X 



1 



17 I J^ -1 



^i) 



a 



Xn 



■^11-^2 -^3/ 
21-^1 -^SJ 



(B.2) 
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The decomposition into K{k), E{k) and n(a, fc) is not unique due to possible 



transformations of the moduh k and a. The ansatz ( A. 7 ) for s = 3 reads 



r(i) 



dxi_ 
dV 



(P) = xj4\p) + «i0^°Hp), i = 1,2,3 



(B.3) 



where ctj = Oj^o- Substituting it into eq. (BT) one obtains the fohowing hnear 
set of equations for the ai 



= Y.^i - 2, 

i=l 
3 

= ^(ai-l)a;2 

3 
= Y.'^i^i ■ 



(B.4) 



j=i 



dx^ 



With the help of eq. ( |B.2D the solution is obtained, which determines the -^rf (p) 
completely 



Ql 



«! 



«! 



1 + 



E{k) 



{xl-xl)K{k) ' 

fi_ A ^ (rrj - xl) E{k) \ 
"" + {xl-xl)K{k)) 



Xn X 



^3 
2 2 



xlEjkY 
xlK{k)^ 



(B.5) 
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